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ABSTRACT 


This work is an attempt at determining the optimum 
Mmeutral axis location in the midship section of a vessel, 
using as a model for the calculations a box girder, subject 
to several different conditions of sagging and hogging 
bending moments, several values of design stress are also 
considered and in the end what happens when lateral hydrosta- 
bic load is considered along with different values assigned 
mo une design stress in tension and in compression is also 


analysed, 
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INTRODUCTION 


It is frequently found in the literature related with 
the strenght of ships such statements as: 
"In general the neutral axis is expected to be at 40% of 
the depth (D) above the base line ... a preliminary estimate of 
the thickness required for the bottom plating (BOT) can be 


obtained from 


O.4D 3 40300 
1.875(——— ) 2240 \ | Ba 
0,6D Lait 


1 S 
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1+ 








( 
950 BOT 


where the Montgomerie's expression is being considered to 
determine the critical stress (SBCRIT psi) and the short edges 


are assumed far apart and simply supported. 








40300 — 
SBCRIT= 
Pee <a 
1 S sS-—> 
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950 BOT _~ a 


a- distance between transverse stiffening members (inches) 
1.875 accounts for a factor of safety 
Tobin's expression was considered in the above formula 


to give the maximum permissible bending stress (LBS, psi) 


B, 
LBs= 2240 \/ AL 





We see thus that this procedure aliows for an estimate of 
the bottom thickness to be obtained after knowing the length and 
assuming a longitudinal frame spacing. What may be questioned 
here is the reason why one assumes the neutral axis 0.4D above 
mae base line and also the reasons behind the choice of a frame 
spacing. 

In the following a procedure to estimate tne scantlings 
of longitudinal members is described; this procedure can also 
be found frequently in the literature dealing with the subject. 

... "For a particular set of principal dimensions , general 
arrangement and frame spacing, plate panel thicknesses for one 
after another of the amidships longitudinally continuous elements 
are selected for an assumed position of the neutral axis at 40% 
of the depth above the keel. The composite midship section assem- 
blage is thentested for the adequacy of its section modulus when 
ranged alongside some standard value, such as given by the 
Load Line Regulations or some function of it. Thus the modulus 
is judged sufficient to limit ship bending stresses to the value 
previously set and incorporated into tne design equations. 

For all but the deepest ships the modulus will undoubtedly 
need to be increased and a small arbitrary increment is added 
to the upper deck or bottom whichever happens still to be farther 
from the true neutral axis location found after the first design 
@ycte,. The shift of neutral axis resulting from the scantling 


gnicrement requires a compensating increment in the opposite flan- 





Ze with the obvious outcome that the neutral axis is shifted in 
reverse with attendant effect on the scantling chosen first; a 
cyclic process is then begun, More trials will be necessary 
before the match of available to required section modulus is 
achieved," 

From this transcript, we may again note the initial 
statement "...for a particular set of principal dimensions, 
general arrangement and frame spacing..." and question whether 
erenot the frame spacing should be chosen as an input or rather 
be left as a variable and be obtained as a result of a structural 
design procedure. It is here aknowledged though that in these 
procedures an acceptable solution can be worked out and each 
new cycle is an improvement over the precedent ones, What may 
be questionable though is whether the first cycle, due to the 
assumptions made at its beginning, has or not a strong influence 
in the outcome leading to a solution among the various possible 
Ones which in the end does not lead to"the optimum design" 
whatever its definition may be, 

It seems that this trial procedure can be looked at as a 
"weakness" of the structural design procedure for the determina- 
tion of scantlings and it is perhaps reasonable to make some effort 
in order to work out a process that could give the scantlings 
mma direct manner and be at the same time in agreement with 
what can be considered "the optimum design", (usually based on 


economic considerations). 
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Nowadays much concern is given to economics in ship design 
and naturally the strength calculations of a ship are also made 
with economical optimization in mind, With this respect, we can 
realize that in ships structures there is a desirability to mi- 
nimize structural weight and some define structural efficiency 
mompeing the least weight structure capable of resisting the 
loads for which they are designed. 

It should be clear that the cost of construction and of 
maintenance and repair are inter related factors together with 
least weight, to be considered when bringing economics into 
Ship structural design, Also other matters such as accessibility, 
weldability and availability of material must be considered as 
well, 

In the following sections an attempt is made to obtain 
without any interactions the scantlings for the plating of a 
box girder. It is hoped that through this simple model some 
helpful conclusions can be made in view of their application 
to a vessel's structure. Each of the different situations ana- 


lized will be exhaustively explained in the following sections, 








Ll 
PRELIMINARY CALCULATIONS 
ime model sed in this section and also tnhrougnout @Whis 


work with only very small additional considerations is a box 


Birger which cross section is pictured in figure 1, 





Pig ae 


As will be seen in the calculations of the following 
Seevions, the dimensional parameters used are of the order of 
magnitude of those considered in ship design. 

This model was chosen due to its simplicity for the 
performance of the following calculations as can be understood 
immediately from its geometrically symmetrical cross section, 

The objective of this section 1s to determine the neces- 
sary formulation to be used in future sections, 

The formulation dealt with in this section is related 
only to the particular model chosen and the geometry of its 


cross section, 
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2(A*d) 
DIBOT= #12 
ZA 
=(A*a) 
DI BOT !'=——— 
=A 


@ A=B*¥TOP+2*TS*D+B* BOT=2*D*TS+B* ( TOP+BOT ) 


DIBOT =Neutral axis location above bottom (in) 


DIBOT'= " " W it (i) 


Leor =Moment of inertia with respect to the bottom 


I = iN il i" T iN u " neutral axis 
3 (A*d): (D/12)*(B*TOP+D*TS) 
I= (a*a*d 5 1~( (A*d))7/( BA) 


Peon ETOP D=TS)/(24)=0SrB8*(TOP+BOT ) ) 
This formula will give DIBOT in inches if D and B are 
expressed in inches; it will give DIBOT in feet in the event 


mio D and B are expressed in feet. 


DITOP=D-DIBOT 
DITOP=D* ( B* BOT+D* TS ) A2*D*TS+B(TOP+BOT) ) 


DITOP in units of inches or feet as explained for DIBOT 
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The following formula will give I in units of (intrt)* 


mor B and D given in inches. 


(D/12)°* ( 3*B¥B¥TOP*BOT+ (D* TS) “+2*B#D#TS* (TOP+BOT) ) 
2 a “a 


3% (2*D*TS+B* (TOP+BOT ) ) 


Next formula gives I in units of (in*ft)* for B and D 
given in feet. 
2 
beD*D* ( 3*B*TOP*BOT+ (D*TS)“+2*B*b*TS*(TOP+BOT) ) 


LS AAA 


2*D*TS+B* ( TOP+BOT ) 


These expressions will be used in most of the following 


sections with B and D expressed in feet, 


NOTE: TS which is defined as the thickness of the 
Side plating is considered throughout this work to be uniform 
and will always be equal to the larger of the thicknesses obtained 


for the top and for the bottom plating. 
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FORMULAS USED IN THE ANALYSIS 


BEAM EQUATION 
Assuming pure bending, the stresses in the cross section 


See cone bOX girder are given by 


BM¥y (/L 


I 






neutral 
surface 


STR= 





x 


Fig. 2 


where y is the distance from the neutral axis to the point 
where the stress is being determined. 
The neutral surface indicated in fig. 2 is composed of all 


ine points which are not subjected to any stress in pure bending. 


BRYAN EQUATION - Plate under inplane load 


a 
(8) *#E*(t)* ——e | SS peg ee 
ea |), | eS 
2a POL) iS) 


When a/S>2 then K=4 
Samem 2£1fen in psi for t andS expressed in inches and E in 


Dei, 
5) =3.1416 
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PLATE UNDER LATERAL LOAD 


S 5 | 
STR3=0. 2 epee \——) a 
(e eeu 





STR3 - tertiary stress 

0 - density of sea water (64 lb/ft?) 

H - head of sea water (ft) 

a - long dimension of plate (in) 

S - short dimension of plate (in) = distance between 
longitudinal stiffeners 


t - plate thickness (in) 


Tne stresses are more severe at points A and B but we are 
particularly interested in the stresses at point B because there 
they will combine with the stresses resulting from bending also 
called primary stresses. 

Mete that for a panel under lateral load, the direction of 
mime Stresses is, at the edges, perpendicular to them, so at B 


they are directed "longitudinally". 
For a/S>2 we have K=0.685 


Replacing all the numbers in the above formula, we obtain 


2 2 
STR3=0,152222*H* (——) 
t 
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The following stress definitions and figures 3 and 4 will 
help understanding the reasons behind the calculations made in 


many of the following sections, 


STRESS DEFINITIONS 

STR1 - SHIP BENDING STRESS or PRIMARY STRESS 

Since this stress is directly proportional to distance from 
the ship's neutral axis, it is substantially constant at any point 
On the cross section of the bottom or top plating when the vessel 
mS in an upright attitude. For distinction, this stress may be 
termed a primary stress or an area stress in this location and 
under such a condition, While the longitudinal direction is imp- 
lied in the foregoing comments, they are equally applicable 


transversely, 


STR2 ~- GIRDER BENDING STRESS or SECONDARY STRESS 

Girder bending stress arises from the reaction of a plating 
stiffener combination to a loading normal to the plating such as 
may be due to water pressure, Following the usual concept, the 
plating acts as one flange of the composite girder and so is 
stressed. Due to shear lag, this secondary stress in the plating 
1S a maximum in way of the stiffener and diminishes with distance 
mom it. 

For the previous case of bottom plating and vessel upright, 
the transverse variation of fore and aft stress will be sensibly 
constant for some distance from the stiffener but will vary 
Significantly through the plating thickness as the plate thick- 


ness is a substantial part of the distance to the neutral axis 
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of the plate stiffener combination, For convenience such a stress 
may be termed a line stress although for a wide spacing of stiffe- 
ners with consequently greater deviation from uniform distribution 
mome inaccuracy Of concept is introduced by using this term, 

For girders continuous over several supports, such as tran- 
sverse bulkheads, maximum bending moments occur at the supports 


and will readily be in the order of twice the field moment. 


STR3 -—- PLATE BENDING STRESS - TERTIARY STRESS 

Proceeding in turn to ever more elemental portions of the 
ship's structure concludes with the simple plate panel supported 
on its four edges, In such a plate under lateral or normal loads, 
bending stresses are created which are variable from point to 
point on the panel and through its thickness with maxima occurring 
at one or the other surface, Since any such stress is effective 
at a point only, it may be called a point stress or tertiary 
eeress, 

Each elemental panel of ship's plating is generally one of 
a repeating pattern, Quite appart from any rotational restraint 
provided by the panel supports consisting of floors, frames, 
Jongitudinals, etc., there is the restraint afforded by the 
continuity of the plating and its loading. For pratical purposes 
then, a condition of edge clamping exists with maximum stresses 
being found at the midlength of the longer panel edges and 
secondary maxima at the midlength of the shorter edges, In both 


cases these stresses are normal to their respective edges, 
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From the previous definitions and figure 3 we may summarize 


Mmitigure 4 the effect of STR1, SsTR2 and STER3. 


STR 
eel 





Fig. 4 


At point A 

STR2 is larger than at point B 

STR2 is acting in the longitudinal direction 

STR3 is larger than at point B 

STR3 is acting in the transverse direction 

STR1 and STR2 are here combined acting both in the 


longitudinal direction 


Ae point B&B 

STR2 may be acting in the tansverse direction 
STR3 is acting in the longitudinal direction 

STR1i and STR3 are here combined acting both in the 


longitudinal direction. 





Usually STR2 is attributed a constant value and can be 
taken into account by automatically decreasing the design stress 


by the value assigned to STR2, 


ae 
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THE BENDING MOMENT EQUATION 


The bending moment equation used is derived from 


BM A*AL/35 ton*ft 


knowing that 


> 


V/35 ton 
and that 
V = BEPHAL*C, 


and using 


we obtain 
OF7> 
BM = ——— *B*T*AL*AL ton*ft 
Sas, 
where AL = is the length in ft 


B = is the beam in ft 
aS ene sarart 32.0% 


Cy = is the block coefficient 


or multiplying by 2240 we obtain 


BM = 1.37143*B*T*AL*AL sHleragiee 


Note: 
V is the displacement volume and {f\ is the corresponding 


displacement weight, 
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MOMENT ACTING IN ONE DIRECTION 


NO RESTRICTION IMPOSED ON THE COMPRESSION STRESS 


Considering the box girder subjected to a bending moment 
acting in hogging BM, we will observe as a consequence that the 
top will be subjected to tension and the bottom to compression, 
These stresses are given by the already introduced formula 


eye 
STR= 





i 


or in this case, being interested on the top and botton, 


BA“DiBTOP 
STH= ——--—____———— 
I 
BM*DI BOT 
SBH= 
I 


iMarmeutis Iloroductory study it was decided not to limit the 
Value of the stress in compression, Only the value of the stress 
in tension, here STH will be requested to be equal to the value 
chosen for the design stress in tension (or smaller). 

SBH will be a consequence of the value of DST and also 
of the S value chosen. 

SO, we will basically have 3 equations 

DST = STH<»TOP Bending equation 


S <> SBH «» BOT Bending equation + Bryan formula 
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OTH 





COMPRESSION ~———|———» TENSION 


Fig. 5 


The three equations that can be obtained will allow one to 
determine the three unknowns which are BOT,SBH,TOP., 

By replacing in the equations already mentioned the expre- 
SSions derived previously for DIBOT, DITOP and J and solving for 
BOT,TOP and SBH we obtain after a lengthy but simple process the 
equations in their final format as we will see. 

Note here that in order to solve the three equations, - 
the thickness of the side plates TS was made equal to the larger 
TOP or BOT; since TOP and BOT were not known and had to be deter- 
mined after deciding on TS, the equations had to be solved twice, 
First the equations are solved assuming that TOP> BOT this 
Simply implying that the neutral axis is going to be chosen alrea- 
dy closer to the top and TS = TOP, Second, the equations are 
solved assuming BOT>TOP thus making TS = BOT. 

As will be seen, only one solution is possible, the other 


being rejected by giving a result that does not agreee with the 
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assumption just mentioned, 


TOP , BOT , S_ (inches) 

Eewe> , DITOP , DIBOT , AL (ft) 
Poi ee scliewecen. & (psi) 

BM (ton*ft) 

I (in*ft)* 


560*BM* (B*BOT+D* TS ) 
DST? STH = 





D* (3* B*B*TOP*BOT+2*B*D*TS* (TOP+BOT)=(D*TS)™~) 


560*BM*( B*¥TOP+D*TS) 
SBH = ————_______.._..._. 
Same denominator 


~ aa 
(s)**(1-(POr)2) 
G = 0. 30396% 


k 


DST Best bD TS 


Making 


we wil 


SBH BewmoetD<S 


9*B*B*D 
A:7-=-—--— 
140*BM 
3*B*D*D 
Co—__— 
70*BM 
jebs* DD 
bE >=... 
140*BM 
1 obtain, by assuming first TOP> BOT 





a¢ 


D*(B+D)*TOP 
DIBOT =——————_———_——__—___—__ 
(2*D+B)*TOP+B*BOoT 
D* (B*BOT+D*TOP) 
Diplo r= 
Same denominator 
(2*B*D+D*D) *TOP*TOP+ ( 2*B*D+3*B*B) *TOP* BOT 
Tl = 4&*)p*p42—2—————— 
Same denominator 
B* (BOT) 32 
ROP a ee _ ae 
(B+D)*G*DST-D* ( BOT ) 
TOP*(12*D-~(C+F)*DST*TOP) 
50) = >» ——_——_————_—__ 


DST* (A+C )*TOP-12*B 


Before going to the equations obtained when BOT > TOP 
is assumed, we call the attention of the reader to the fact that 
the values of BOT and TOP have to be obtained from the two above 
equations, Substitution of one into the other would result in 
an equation of degree higher than at least 5th degree, unsuitable 
mor airece solution, As will be later exemplified, the solution 
is here obtained by assigning values to either one of BOT and TOP 
and through the two equations, we obtain the respective values 
for TOP and BOT, The ploting of these values BOT vs TOP will 
mye one line for each of the two equations. The intersection of 
the two lines thus obtained allows the solution to be read from 
the plot. 

This solution will be acceptable if the intersection lies 
in the sector TOP> BOT otherwise the intersection must be 
me jected. 


Assuming now BOT> TOP we make TS = BOT and obtain: 


r 
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D* (B¥TOP+D*BOT) 
PL — 
(2*D+B)*BOT+B#L0OP 
D*(B+D)*BOT 
DITOP = 
Same denominator 
(2#*B#D+D*D) *( BOT) ©+(2*B*D+3*B*B) *POP* BOT 
I = 4*D*pD* 
Same denominator 
(B+D)*( BOT) 3-p#*G*DST*BOT 
TOP = 
B*G*DST 
12*(B+D)-(A+C)*DST#*TOP 
BOT = 


GEE a | eo SEL 


In this form the equations can be used in a very simple 
computer program that will give the set of values TOP vs BOT 
Suitable for ploting. At this stage we are hoping that if the 
intersection obtained in the situation where we assumed TOP> BOT 
is acceptable,the two lines resulting from the case BOTYTOP will 
meu intersect or if they do that the intersection must be rejected 
@ue to the fact that they will violate the assumption from which 
they were derived 

As can be observed from the several cases run in the compu- 
ter, the case where two acceptable solutions exist never happe- 
ned and it can only happen if TOP = BOT where then both solutions 
are acceptable and coincide as in the case of one of the follo- 
wing runs, 

In the following pages the computer printouts are included, 
The first run is explained in detail to make the interested rea- 
der familiar with the problems encountered in determining the 


intersections and the respective thicknesses. 





ay 


The other runs are here presented as well to serve as 
proof that a simple solution can always be obtained and also to 
meme am 10ee@ of the range of thicknesses expected in each case. 

In ali runs only a bending moment acting in hogging was 
considered. 

The program used is a direct application of the eauations 
here derived to a particular set of values chosen as the input. 
From the results obtained, some conclusions were drawn that a- 
liowed for the improvements made in following sections, 

No further explanation on the program shall here be given 
because it is self explanatory in the comment statements 


included. 
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From these lists of thickness values printed by the com- 











maver, we can see that the acceptable ones are those marked 
with a Siem, Frlotting these corresponding pairs of values 
meoperly, as is next exemplified, allows for one to read from 
mie intersection of the two lines, which falls within the range 
set by the assumption that originated it, the values of thick- 
nesses at the top and at the bottom which satisfy the model 


we Started with, 
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We can see that the thickness values are easily read from 
the intersection obtained and they are acceptable because they 
fall within the range assumed to derive the formulas that ori- 
ginated them namely TOP> BOY. 

Next will be plotted the case BOTD>TOP and as can be seen, 
peain an intersection 1S obtained which in this case is not 
mecepvablie because it occurs within the range TOP> BOT contrary 
to the assumed BOT> TOP that originated it. 

With a little more experience, this can be seen immediatly 
from the computer printouts, thus avoiding loss of time doing 


splot that will lead to a nonacceptable result. 
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Now having obtained the values for the thicknesses from 
tne plots, we verify them through the use of the next program 
Seased on the formulation previously derived, As can be noted, 
we obtain an acceptable value for the stress at the top (since 
we are considering a hogging bending moment) which we said 
before in the derivation of the equations should not exceed the 
@esien stress in tension, In this case the design stress was 
DST=20000 psi and we see that the stress at the top STH=20037.5 
the 37.5 coming possibly from inaccuracy when reading the thi- 
Gknesses from the plot, and computer roundoffs. 

Note that the stress at the bottom, since there was no 
limit placed on it except the Bryan formula to resist buckling 
meused by compressive inplane loading, is attaining a too large 
Value which iS unacceptable and unrealistic, Here again there 
meea be no difference between SBH and SBCRIT and the difference 
that can be noted is also due to inaccuracy in reading the thi- 


cknesses from the plot and computer round-offs. 
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SUMMARY FOR THE PRECEDING CASE 


Length = 500 ft 

Beam =AL/ 8 

Draft = B/3 

Depth = ALAS, 

DST = 20000 psi 

5 = 10 inches 
Material -~ STEEL 

TOP = 0,468 inches 
BOT = 0,166 inches 
Sy ive! = 20037.45 

Sleval = 30434, 89 
SBChET = 29886, 80 
DIELOP = 22,06 ft 
DIBOT =o Ow 


We can already notice here that the stress in the bottom 
is attaining a too large value due to the fact that, since no 
limit was imposed on it, it is dependent only on the chosen 
Maiue for S and the resvective capacity to resist buckling pro- 
vided by S and BOT while at the same time stressing the top plate 


to its maximum allowed DST. 


In the following two cases only the spacing S between the 
longitudinal stiffening members will be allowed to change and 
the results obtained are simply stated as previously exemplified 
and were determined in a similar way. 


From the results obtained we can note the effect caused 


by the change of S. 
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YeOos co} 37 ag° 2s ZeqgkZZ 
Toe dew 10.4925 Suc Mee 1°@7ee 
14? gic Ae1erD BaT4s P2996] 
TAD YAS PG A18KCY RATh&e 729731 
TnPu= G.Prr? 2074s Feri et 
ToCas 722574 5C 74> 7°7170 
bag one 7 eo ZPPRD Robe 7 *h6P4O 
TayPas 003572 Bogs = Zer aig 
(“Pac eh AGT ROT4Ss 7°4379 
TOPE 2 eh 502 R074 FeBnag 
TOP,& & «eS AGie BTSs 2°2519 
Teas Gea Sae QOT4&s YAM 001 S) 
TAPAS FG ebFDA Qnvras 720 45k 
TOPPA Geh570 BOT4= #09272 
TOPas GC e7AL? ROT4&e #704 2%e 
Wier ¢ © 7584 = 6 ie Ae eh TB) 
quel. FeFARD RO Y4se mee 3763 
TOF AS 705592 BQT4= #7+3993 
TeSys 7 eGAPD ROT4= "7049203 
7 Oeae= 7 eP509 BOTS4= #2 eBrh4 
TOP AS 172279 BO7Ss, ~7+4784 


END 


VIC OPERATING SYSTEM VERSION 1 REVISION 012 @370@4/74 GENERATED @5/06/74 7? 


JOB -ANOLING CHARGE > .,35 / JOB 35 
73 LINES PRINTED PRP 2 1.°5 7 K LN ©2? 
67 CARDS READ me? 7 K CD °1¢ 
me “LOTTER VECTCRS Se wet 6/1 RC °77 
Mee PODEL 72 SECONDS P5.407 7 WCUR °14 
Mme ODEL S2 SECGN®ES $12.87 / #OUR 7? 

TOTAL CHASSE $ *Ri 


Memeo 1273; LOGGED CUT 75706/74 15:24654 8 13429 LEFT AFTER 18 EOQGINS. 
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STRENGHT Of SHIPS©ON THE NEUTRAL AXIS LOCATION PAGE 
W204 C ~A_* IS THe LENGHT 
PPA AL=5A06 
pr ac AREAL SR 
WIR C sORA® 1S TE DRAFT 
-18 DRAZR/3 6 
go 2% NEAL/9-° 
"39 Ee «POl* IS THE POISSON'S RATIO OF THE MATERTAL 

22 C FOR STEEL 

22 PO] she 3 

°33 C «cm IS THE YOUNG'S MODULUS OF THE MATERIAL 

zi a FOR STEEL 

Ne Sa = 1 eG 

We 4 ~L$T* IS T-f DESIGN STRESS 

ig Sule > (* Soils 

ocal - eSe IS Tat SPACING OF LONGITUDINAL STIFFENING MEMBERS 

= SF3o 6 

se ‘a eteeu TS THE RENDING MOMENT 
1950 WV (RRNA AL AL He FE) /( 3250 #350) 

a ge SE 39S=S5FSe (1 --PClmP D1 YE 

Fe C eTOP= TS THE TMICKNSSS OF THE TOP PLATING 
Sok c e207 TS Tek THICHNESS OF THE BOTTOM PLATING 

— TPE e 4425 

43¢ NOTH 6412 

a SSSA Meee MM (QR RAGT He DHA TURP ) 

— “ENED (SRI KE HK TOOARDT AC eee Du TOP ( TOP+HOT) 4(040% TO? «TOP ) ) 
140 c «Tle IS T-5 STRESS OX TOP OUE TO HOGGING 
PLGe STYsA/OAN 

15a Bret. 6h ySTt 

(7h ee ORMAT( ' ' Pex, "STs! F117 +3) 

“5c So er. oer eg] OP + AT OO ) 

:= C eo Soe SEK ESS ON BOT OUE TO HOGGING 

oa 3M t= OZ ODEN 

cA RT T&1S,9215aH 

2 DM TCR AT(* ' 2x, 'Gate (51043) 

na - ‘a Seer tA 45 TA Gar ICAL STRESS ON SOT 

eg eee Il =SC iment 75 

Date RITE(R, {Lo )SUCRIT 

225 Sok me ee! FC GLK "SBA CRITSEIE Ce 3) 
Pees 2 S1De FP TOP + Ral TOS4 9 OT ) 
Rohe c -BETOP. TS THE DISTANCE OF THE TOP FROM THE NEUTRAL AXIS 
P2aR NT TCOaDe (he RAT + De TOAD) /P 
Poor ae’ mel > Oe) DTT OP 
ro 4a eee ie eAT(* 8 4X, POT TOPS F146 & | 
Poa C Om ee TS Tee DISTANCE OF THE BOTTOM FROM THE NEUTRAL AXIS 
pecs NDIBNT=HED#(KReTNF4+De TOO) /P 

wae WRITE(Ss5V)O0160T 
230% ac) comeicat' , Ww, 'CWGOTEIF 1A. 4) 
C 76g, END 
PetS) ev 6 324CV) AL g32eC lv) B G334(V} DRA W33CCV) D 
(yj) Por +34C(V) £ €GV7(S) «R CSG ty} CST O264 CV) o> 
te (CV) bs W37C(V) G G3Re CV) TOP A39@(V) BOT 2394 (VJ A 
Pstv? cen M3B4(V) ST. ?174(L) 67 O772CS) @1 23330V) C 
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ReReNGh! Of SRIPS*®ON THE NEUTS AL AXIS LOCATION PAGE Cc 
|ctv? S&b CLEA2TL) 92 Q2C. CV) SBCRIT weeaetl) 110 e3C4CV) P 
|e ve err oe CCAE(CL) 4¢ CSCC CV) Sb 1E0T Voea(lL> “50 QHAAATS) «-V 
rEQ e 
ci 
OGmewr LABELS. 
Pay eMATN® CAAA oV eCfO@ «COMP CAGE @] C447 oR C938 ee 
DhAn £6 CLR eA CORK eo“ES PRX3e2 eW C63K BP e7CO Als 
1253 ALCG CRFR eRARG COHK 6S PA3A AEXP 295C e ERENT COCA HER 
295 +0 CASe eV SPRA 
TRY" FCINTS: 
bet’ CHE oA C52P ALOG PABA AEXP 263K EXP C70C@ AIN 
jes2° °% eroa¢ .CoMn CERA $46 PRFBR eRARG COIOLA FB C954 65 
950 «ZERO ©C95C -ERC.T C9OSF 20 2984 «MES CABE VU CAKH oV 
Pee ¢! 
WHOMFPLOCKS: 
Bee 
ia 
! 
PEF INED SUBROIJSTINES: 
wo 
i} 
| 
pSFER mene' SS e776 
Sues ROGTS: 
| STH= 2°C37-113 
SBH= 2°737+-113 
SRERIT= e€° 575.346 
ol1TerP= C7e7ITSE 
| OIBOT= 27+7778 
0 
VIC OFPERATI?NG SYSTEM yYERSION 4 REVISION Ale 03/04/74 GENERATED 25/06/74 V 
JOB HANOLING CHARGE $$ .35 / JOB °35 
SS LINES PRINTED PRO © 1.25 / & LN oe | 
5& CARDS REA ¥1.oy 7 kK CH 0VR 
Me PLOTTER victors 25 / 17 Ae ©? 
BeerecUlt 7%} SECONDS $25.77 / HOUR 2793 
MeeroOct §&4 sECCNOS 212.57 / KOYR al) 
TOTAL CHAwEe = Ta 8 


Mies $4731 LOGGED OUT 





D/GA/74 13.466 


S diel, LEFT AFTER 21 LOGINS. 
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SUMMARY 
Length =Ee OOM ing 
Beam SG 
Draft = 18 
Depth = AL/9 
Si = 20000 psi 
3 =. 0! 
Material - STEEL 
TOP =O AS 
BOT =O bi 5! 
Saul = 20037.11 
SBH = 20037.11 
SBCRIT = 20505.43 
DELOr = a ae 
DLSOT Smee (ao. 


Note that this is the case where two acceptable solutions 
exist although they end up being the same because TOP=BOT, 

Note that the neutral axis is located half way between 
top and bottom and the stress is the same at the top and at the 


bottom, as it was expected beforehand, 
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STRENGHT OF SKIPS*ON THE NEUTRAL AXIS 


iA 


a0 
a1 


4? 


‘in 
2 


E60 


aul 


70 


LOCATION 


*eALe JS 

AL #5@He 

B=AL/&e 

*DRA* JS 
NRA=R/3 6 
D=AL/G e 

*POT*e IS 
Sr STEEL 

POT SQ«3 

WEES THE YOUNG! SeeMIgDULWG OF THIBMMATERT AL 
Ree Teck 


wide OG 


THE LENGHT 
THE DRAFT 


THE POISSON'S RATIO OF THE MATERIAL 


CSQee 

w-STm IS THe OESIGN STRESS 

JST rev Cae 

*Se TS THE SPACING GF LONGITUDINAL STIFFENING MEMBERS 
S=i2t. 

save IS THE BENDING MOMENT 


wMe( Re DK AZAL a AL ©4275) /(35e"35 6 ) 
ae deere ee tee OT ero] ) /E 
A>Germe ene Ds (447 en AM) 
raQee tude Ds (7202 A) 
Se im a me 
Berk SS 8 TCR 8S5SUmra LARGER THAN SOTTOM TRICKVESS 


',eSX,' TOP ASSUMED LARGER THAN SOTTOM' ) 
eTOPs ITS Tuk THICKNESS OF THE TOP PLATING 
e7CMm IS 726 THICKNESS OF THE BOTTOM PLATING 
SOT tals 
~CO*SReAQT, RROTL ROTI ( (24D) #GeDST9 Oe BDTL*ROTL ) 
Maer CS remy TOP. AOT + 
Famcere ("8 PSK, I TOR SIE Ge hy OX, MOTI ='h 9.4) 


wey Ren] fw ig AS 
BOO lista ee GO. To 3c 


TG ean 
TOE P= 2? 
ALT P= OPP (1 Pew NH (THF) KR OST*TOP2)/S(DST*(A4+C) *TOP2"12+%*B8) 
WReTTEMSs,aPyvTOPPrs,a9T2 
BORMATEt §',PaK,'T9IOP=tIFO. Nyse os Ue tenpae 'F9.4) 
FUP P STOP DO HG ee AP 
me sOP Pei 2.1. )G0 TD 52 


oe « 3: 

Oe ie inte LS 

TY ICYNESS CF BOTTOM ASSUMED LARGER THAN 
wey 5 (5, 4G) 

CORMAT('L',25X,'ROTTOM ASSLMED LARGER THAN TOP!) 
SOT 3=1« 

FOER=A( L4H) *AQOT3*e ROT Ae BOT 3-D¥eGeDOST*801T3)/(B*G*DST ) 
wWRYTEI(5,77)TOPZ»,RQATZ 

FORMAT (! ",25Xe'TOPZStF9e4,5X,'BOTI='F 944) 
BOTASKOTA+G 675 

TF (BOCT3*+GTe2e) G0 TO 8A 

Ge TO 61 


noe 
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STRENGHT Or SHIPS=CN THE NEUTRAL AXIS LOCATION PAGE C 

4Ae 8A TOPSSt ee 

Noga R1 nOT4#2#(12e¢*(B+0)=(A+C) *¥DOSTRTOPG)/( (CHF) &DST) 

Pu4e WRITE (5s 9@C )TOP4S*POTS 

13472 9G FORMAT(!' ',2P&X,!'TCPyactFGeh4s S5SX,'BOT4E'F 9.4) 

u9A TOP4Se& TOP 444.272 

24 AG PPC UCP Ae Gla. ) GO TG 7aA¢ 

7L GER GC FTG ei 

24BC 420 CONTINUE 

ALAC Ep 

2@(S) «VU S4CH0V) Al. C4CR(V) B Q4NUWCV) DRA C4D8CV) DO 
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pO) 12 PYPY CS) a] C53CCv) BOT1 0176CLI 11 ZO4erv 9: 7 Oka 
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Vee «=2TOP3 Poer iL. 76 AKEPPEL) 82 @574 CV21 TORS A4AACLI 81 

7&Cv) BOT4 me7ett. 9f C@PCttL2 122 CY@ALS) «V 

| x£0 L 

1 

SOGRAM LABELS: 

| 2076) *MAIr~ 2C*CE eV 2@43C «COMP PC1A @I PSEC «R 2BiN4 «Zk 
crhbe 4G C6°C eh CREG «MES eCSDE ew c7EA EXP e96C Al 
c6D0& ALQG ZAAG eRARG CALE 0&5 27E6 AEXP PROB eERCNT CACH $8 
CBA 00 2RBDE eV 3N364 

NTRY=PQOINTS: 

| CECE el eC4Pe ALOG P7TEG AEXF C7EG EXP e964C AT? 
29DE ew Pa3C «COMP CAEP $6 PAAR eRARG ZADE $8 CRAG 5 

| cPG@S& -ZERO CARGR eERCAT CRPA 0 2B60 eMFS PREA eU 2C1if eV 
Ge =! 


INDEFINED SUBROYTINES: 
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RENSFER ADDRESS 2076 
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Note that in this case by simple inspection of the 
coordinates of the points to be ploted we can see that the 
mntversection will occur outside the acceptable range, 

Note also that in order to obtain the coordinates 
of the points within a range of interest where the intersection 
May occur, different ranges and increments are considered in 
each program, This could be improved if at that time these 
printouts were available and a better idea of the thicknesses 


that should be expected existed, 
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SUMMARY 

Length =m) 
Beam = Al/ 
Draft = B/3 
Depth =7AL/9 

DST = 20000 psi 
S = 120" 
Material ~ STEEL 
TOP = 0,03" 
BOT = 15 
Sfa = 19771. 88 
SBH =e iO 
SBCRIT = 9788.37 
DLLOP = 37.42! 
Dior = 13713! 


Note in this case the 


large value chosen for S and 


associated with it the large value required for the thickness 


at the bottom to avoid buckling which nevertheless does not re-~ 


Sist stresses larger than 9788 psi in compression, 
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If we look now at these three cases, all for the same 
material and set of dimensions (same bending moment) we can 
see how the thicknesses and neutral axis position changed as 


the value of S was changed. 


BOT SEH SBCRIT 
Miehnes| inches jinches psi psi fate 
“ | ie 
: | q 
4 | : 


Ineschis table, we can more readily make conclusions on 














the influence of the value chosen for S on the neutral axis 
location, and also a quick look at the resulting values for 
the thicknesses may give an idea of the weight associated with 


each case, 


In the following cases the dimension parameters are 
changed (increased) and consequently the bending moment and the 


thicknesses of plating will be higher. 
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STRENGHI OF SHIPS*ON THE NEUTRAL AXIS LOCATION PAGE iL 


Cc *Acw IS THE LENGHT 
ALSlAh. 
BeAL/5e75 


C *ORA« JS The DRAFT 

SRAZB/3¢3 

D=AL/14% « 
a sPOI* IS Tee POISSON'S RATIO OF THE MATERIAL 
a BOP STEEL 


POTHO3 


a foe > Wee YOUNG'S WVPODULUS OF THE MATERTAL 
a FOR STEEL 
2ST eH 1 BM Oh 
a *>ST*« TS Tower CESIGN STRESS 
DHT = LP stated © 
a SSR 15 THE SPACING OF LONGITUDINAL STIFFENING MEMBERS 
ZeAue 
3 meMS (5 The 3ENDTSG MOMENT 
2] ~ 
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ee bales (7d enh) 

ei <i BD, ( 1 aes & ae) 

PES, me 

10 FORMAT('1!,25%,' TO" ASSUMED LARGER THAN Be OMe) 
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52 CONTINUS 
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WRITE (5,40) 
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Length 
Beam 
Draft 
Depth 
DSP 

S 
Material 
TOP 
BOT 
STH 
SBH 
Seca it 
DITOP 


DIBOT 


Ve 


SUMMARY 


1000 ft 
AL/5.75 
BY) 348 
AL/14 
18000 psi 
BOE 
STEEL 
baal 
Omogt 
MICS LOZ 
a re 
46323. 82 
OR oS 
pie 7 


We might say here that the value chosen for S is small 


and this accounts for the compression at the bottom to reach 


Such a high value, 


Peeune mexeecace, the value chosen for 5 is increased to 


60 inches. 
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SUMMARY 
Length = 1000! 
Beam =e iNley sey S 
Draft = B/E 
Depth = AG/ 14 
DST = 10000 psi 
S = 60" 
Material < (SHO eM E, 
TOP =the Ot 
BOT = 1 1 
STH = 18017. 24 
SBH = 39055.43 
SBCRIT = 39843, 44 
DITOP — 2 Go 


DIBOT = 48, 88! 
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Sinee sometimes it is desirable to use aluminum alloys 
mo the structure of ships, due to their light weight, high 
corrosion resistance and non-magnetic characteristics, in the 
BoLlowing runs, the material considered is aluminum, The 
necessary change is made in the input to account for the pro- 
Derties of this material. 

When considering aluminum as the structural material, 
we must Keep in mind that the aluminum hull will be more flexi- 
ble than a corresponding steel ship. The modulus of elasticity 
of aluminum is about one third that of steel, and if scantlings 
were identical, the deflection of an aluminum ship, under a gi- 
ven load, would be three times those of a steel ship. 

However, to obtain a factor of safety on ultimate strength 
mae! to that in steel, the area of an aluminum member in direct 
meress, or the section modulus of an aluminum member in bending 
would be about 1.5 times that of the corresponding steel member 
(based on an ultimate strength of about 60000 psi in medium steel 
and of about 40000 psi across the welds of a welded aluminum 
eeructure ), : 

A section modulus of 1.5 times that in steel, with the 
same depth being considered in both cases, would result ina 
deflection in aluminum about twice that in steel, However, the 
lower modulus of elasticity of aluminum reduces its ability to 
resist buckling , and for the same factor of safety both the 


Section modulus and the moment of inertia in aluminum would in 
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general be somewhat more than 1.5 times those in steel, so 

that the deflection of the aluminum structure would in general 

be somewhat less than twice that of the steel structure. 
Wieweeealverwmac! leculen would reduce Che natural frequency 

of hull vibration, and sowould affect the interaction of slamming 

stress and bending moment stress, It could also affect shaft 

bearing pressures. It has been proposed as an arbitrary basis 

Of design that aluminum ships be designed for a deflection about 

1.5 times that of steel ships. This would require a moment of 

Baertia of twice the value of that in the steel hull girder, 

and would automatically result in a section modulus substantially 

more than 1.5 times that in steel, and a factor of safety cor- 


respondingly greater than in steel, 
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SUMMARY 
= 1000! 
= AL/5.75 
=—wBy 3a3 
= AL/14 
= 18000 psi 
= 120" 
- ALUMINUM 
= 4.27" 
= 2.98" 
er 9 DO 
= 22852.76 
= 22768.12 
= 31.43! 
= 39.99! or 40.00! 
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Looking at the following tahl€ we can also observe for the 
aluminum the effect of changing S on the neutral axis location 


and the compressive stresses resulting on thebottom plating. 


CASE S Ler BOT SBH SBCRIT PEO e BAS OMe 
inches | inches] inches psi psi Ct ag 
x 
10 oa Ons63 2000 49000 Poo oem 


eS 


4 120 4.27 2.98 23000 40,00 


(*) This difference may result from inaccuracy in reading 


the plot, 


Here again the values of TOP and BOT may allow for a quick 


evaluation of the weight relationship between the four cases, 


One run was also made in the case of the aluminum for 


AL=500 ft and the results are here included. 
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SUMMARY 
Length mo Or 
Beam AL/8 
braft B/3 
Depth AL/9 
DST 18000 psi 
S 30" 
Material ALUMINUM 
EOP Di ie) 
BOT GAce 
STH 18089, 66 
SBH 1457814 
SBCRIT 14767,85 
DITOP 30.76" 
DIBOT 24,79" 





HOGGING BENDING MOMENT=SAGGING BENDING MOMENT 


DESIGN STRESS IN TENSION=DESIGN STRESS IN COMPRESSION 


From the preceding section we can understand the need to 
But some limit on the value of the compressive stresses, and per- 
haps also the inefficiency of considering the value of S (the 
spacing between longitudinal stiffening members) as an input to 
the problem, 

In this section we will consider a bending moment of the 
same magnitude in both hogging and in sagging and the stress 
in either tension or compression will have as upper limit the 
gdesign stress, 

Since the bending moment is the same, it is intuitive that 
the thickness required to resist tension on top and compression 
on bottom (buckling) when the hogging moment is acting will 
also be respectively the thicknesses to resist tension on bottom 
and buckling on top when we consider the sagging moment, We 
mould end up then with a thickness on top required by the 
hogging moment TOPH, a thickness on bottom required by the hog- 
ging moment, as well as a thickness on top required by the sagging 
moment TOPS and a thickness on bottom required by the sagging 
moment BOTS, 

Since HM = SM we must obtain 


TOPS = BOTH 


ti and BOTS = TOPH = t2 
eoeWemwt@et Cnd Up considering the largest tl or t2, to 
be the thickness of the top and bottom simultaneously. 


Assume for instance that t1>t2 then if we decide on tilfor 








2 i 


thickness of the top as well as at the bottom, we will 
verify that we have overdesigned because TOPS = t1 
Corresponds eOnsol.s — Ccowaten plugged into the expression 
for I gives a smaller value than using TOPS = t1 with 

BOT = tl tnat will give a larger value for I; this larger 


I value when used in 
BM*¥Y 





SIR. = 
it 


will give a smaller STR value than the design stress and 

we could use a smaller thickness, may be between tl and 

t2 for both the top and bottom and satisfy the DST requi- 

meme iy. 
This means that we must incorporate in this analysis the fact 
already known that when BM (hogging) is equal to BM (sagging) 
we will have to consider the neutral axis half way between top 
and bottom which means that the thickness of the top shall be 
maual to the thickness of the bottom for a geometrically 
Symmetrical cross section, 

Let's determine the formulas that will best serve this 
Base then: 

HM = SM = BM 


DITOP = DIBOT D*(B*t+D*t )/(2*D*t+B*(2t ) ) 


t = TOP = BOT = TS 


3*B*B+D*D+4*B*)D 


tH 
il 


C*DeDaie 





We aye) 


BM*DITOP BM*DIBOT 
Now we can write STR = ————- = ———— = DST 
I I 











Ve 


and the instability formula 





1-POI*POI 
oka 0, 30396* ———--—__——. * DST 
E 


Where S is the spacing between longitudinals. 

We see thus that we cannot have more than two equations 
and so we must seek for two variables, A quick inspection and 
One can see that knowing BM, B, D, the material properties and 
the design stress DST we have left to determine what should be 
the thickness of the plating (t) and the corresponding frame 
spacing (S). Solving the equations with this in mind, we obtain 

BM = 1,37143*B*T*AL*AL lb*ft 


ley uland 


and from Sui = we obtain 





a 


BM*(D+B) 
DST = 





LeD*t (3*B*B+D*D+4*B*D ) 


which gives 


BM* ( D+B) 
t =———— 


4L*D*DST* (3*B*B+D*D+4*B*D ) 


and also 





0. 30396*(1-POI*POI)*DST 
E 
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These formulas will give for a prescribed set of dimensions 
(ft) and a particular material subject to bending BM (lb*ft) the 
thickness t (inches) and the spacing of longitudinal stiffening 


members S (inches) required. 
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HOGGING MOMENT AND SAGGING MOMENT OF 
DIFFERENT MAGNITUDES 
THE DESIGN STRESS IS THE SAME IN TENSION AND 


IN COMPRESSION 


Let's consider now the case where the magnitude of HM is 
marterent than the value for SM, 


We have in this case: 


D* (B*BOT+D*TS ) 


IMO yp 
2*D*TS+B* (TOP+BOT) 
D* (B¥TOP=D*Ts) 
Se. as 
SAME DENOMINATOR 
3%B*B¥ TOP*BOT+ (D*TS) <+2*B*D*TS* (TOP+BOT ) 
I = 4*p*p* 
SAME DENOMINATOR 
HM*DITOP SM*DITOP 
a BTN ae 
i i 
HM*DIBOT 
= SM*DIBOT 
I Oo elas 
i 


BOT = \c \\sBE Or \ SES 


s<*(1-POL*POT) 
where GS 0 oe Sige 





iD 





oe 


Looking at these equations, we know that the result obtained 
for the thicknesses is going to determine the value of DITOP, 
DIBOT and I. If furthermore we restrain the stresses from being 
larger than DST we will undoubtedly want them to end up being 
DST otherwise if they are smaller we are not using the material 
efficiently. 

We thus have 
WOE, & 


Dor = STH fall eel Dole =— Sls = she Kkt Where ai 


tt 
tt 


| 


DST = SBH = HM*Ke DSL =opoo= ol~Kz where K2 DIBOT/I 
This tells us that if for instance HM> SM 
we should be considering DST = STH = HM*Ki 


and DST Dede We Ke 


because the other two equations will determine lower stresses 

and are thus within the requirement of not exceeding the 

design stress, being automatically satisfied. 
mow from these equations we can see that if the design stress 
required for the top, (in tension for the case of HM>SM; note 
that for the case SM>HM the other two equations should be con- 
Sidered and the conclusions are similar), is going to be the 
Same one as required for the bottom ( in compression ), we must 
have Ki = K2 and so DITOP = DIBOT and the neutral axis wil? 
again end up being located halfway between top and bottom for. 
&@ geometrically symmetrical cross section, 

We can here use the same two equations derived in the 
previous section to determine (t) and (S), The moment to be 


considered here will be the larger of HM or SM instead of BM, 


i 
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HOGGING MOMENT LARGER THAN SAGGING MOMENT 


DESIGN STRESS IN COMPRESSION LARGER THAN DESIGN STRESS IN TENSION 


Thesituation we are analysing in this section 
HM > SM 
DSC > DST 
is perhaps more interesting than the previous ones because its 
study provided us with a much simpler method for determining the 
thicknesses and the neutral axis location; we call here the atten- 
tion of the reader to fully understand the steps of this section 
because the following sections till the end of this work will use 
this same method of analysis applied totheir own particular case, 
Indicated in the next diagram are the ranges within which we 
May stress the material in order to satisfy the requirement of 


mot exceeding DSC or DST. 


at———— SC DSL 


Oo 


COMPRESSION TENSION 
§———$ 


Let's consider now that we will try to determine the best 


position for the neutral axis by first gaining some insight into 





SA 


meee comepry Ol tne Stress GQi%ecributions and the stress limits 


already defined, Unacceptable 
| | t 1 
~—— DSc -— DST» cc 
| : z 
: 4 
a 
DITOP | va 
Ca | 
NEUTRAL AXIS 
i —-—-—- — = ene 
| 





Acteptable Sia | 


From the previous diagram we may easily understand that 
if we decide to stress the material to its maximum in tension 
and compression when under the effect of a sagging bending 
moment (SBS - STS) we end up with a neutral axis closer to the 
bottom due to the fact that DST < DSC and we see that when consi- 
dering the hogging moment the only possible alternative left 
is to limit the stress caused by this moment in tension on top 
STH=DST; in this case the compressive stress on the bottom SBH 
would automatically fall inside the acceptable range, 

If for the hogging moment we choose to allow the stress 
in compression on the bottom to attain its maximum value SBH=DSC 
we see from thediagram that we are forcing the STH to fall out- 
Side of the acceptable range. We are left in this situation with 


an acceptable solution so far and we.must investigate it further. 
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At this stage we have 


IS 


: 


SBS=DST 


We may see from the diagram that we are asking for SBS=STH= 
=DST to have the same value and at the same time to have the 
meutral axis located closer to the bottom DITOFP> DIBOT, 


We note that the equations that give STH and SBS are 


STH HM*DITOP/I. 


and S188 SM*DIBOT/I 


and since I is the same STH can never be of the same magnitude 


as SBS because 


Sth = HMpitOr/i > SBS = SM*DIBOT/1I 


both because DITOP is being set larger than DIBOT and also 
because we are considering HM>SM to begin with, Thus we are 
left with the alternative of making DIBOT larger than DITOP 
to make it possible to satisfy 


Stee —-eEM=pinOrP7i = SBS= SM*DIBOT/1 
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Thesituation would then be 


DBC SLHEDST 

, 

| 
ee hae =e 

| =i 

| , 

) | 

| , 

| | 
DSC=SBH SBS=DsST 


Besides this is only natural because what we are doing here 
is to make the material to be stressed to its allowable limit 
while resisting the larger of the two moments; hogging in this 
case, This situation will be worked out here as an example and 
it seems to be quite clear that similarly, cases where SM>HM or 
Mol >DSC or other situation like DST >DSC with HM>SM or 
DST< DSC with SM>HM could be solved applying a similar reasoning, 
50, having HM>SM and DSC >DST we should consider the smaller 
design stress,in this case DST,and have the material stressed 
®©O this upper limit. Why?,because if the larger is considered, 
DSC in this case, there is no way we can have the material stre- 
ssed to this higher limit in both hogging and sagging while con- 
Sidering the smaller we may attempt to do so, This can be better 


understood from the two next diagrams, 
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DSGeors DSl=S cH 
| | 
. | 
| i 
ee ee eee 
: | 
| 
| 
| 
| ’ 
. | Too large 
| 1 SBS 
DSCFSBH 
DSC=STS SME Toolarge 
, , TH 
| | 
| 
| 
| | 
! | ONA, 
| 
| | 
| | 
) 
WS GS S15) a! DST=SBS 


So, considering the lower design stress, in this case DST, 
we will stress the material to this maximum in both hogging and 


Sagging thus giving 


STH = HM*DITOP/I 


DST 


SBS SM*DIBOT/I DST 


Note here that we already know that the neutral axis will 
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be closer to the top DITOP< DIBOT which means that we must obtain 
a thickness for the plating at the top larger than at the bottom and 


so TS=TOP which means that we will be using the following for- 


mulas 

D* (B¥*BOT+D* TOP ) 

DITOP = —————--—___—_— 
P 

D*TO P* (B+D) 

DIBOT = 
P 
ux D“*TOP* ( (3*B24+2*B*D) *BOT+(D*+2*B*D) *TOP) 

LS ee Ss eee 


le 


where P =(2*D+B)*TOP+B*BOT 


we have also 


STH = DST = HM*DITOP/I 
SBS = DST = SM*DIBOT/I 
SBH = DSC = HM*DIBOT/I 


From these three equations we may combine the two first to obtain 
HM*DITOP = SM*DIBOT 
and combining the second and third we obtain 
DST*HM = DSC*SM 
After combining these two equations we obtain 
DST*DIBOT = DSC*DITOP 
This equation could have been obtained more easily just 
by looking at the diagram we are considering and notice the 


Similar triangles 
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IDNs 





DSC 


Note also that what we obtain from this equation is to 
determine the N.A, location by making DITOP and DIBOT proportio- 
nal to DST And DSC. 


Proceeding with this equation we obtain: 


DST* (D*TOP*(B+D))/P + DSC*(D*B*BOT+D<*TOP)/P 
TOP*(DST*B+DST*D-~DSC*D) = DSC*B*BOT 

and making 

AC = DST/DSC+D* (DST-DSC)/(DSC*B) 

we obtain 


BOT = AC*TOP 


me need now an equation to stress the material to its design 
memit and this can be done by either using 


STH = DST HM*DITOP/I 


or SBH = DSC = HM*DIBOT/I 
meom these two equations I prefer to use the secondespecially 
because it will simplify the calculations that follow as can 


be understood from the simpler expression that gives DIBOT as 
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compared with the one that gives DITOP. 
We will then consider here 
SBH = DSC = HM*DIBOT/I 


which gives 
HM*D*TOP*(B+D) 
DSC). ————— 
hep?sTOP( (3*B°+2*B*D) *BOT+(D2+2*B*D) *TOP) 


Simplifying we obtain 


HM*(B+D) 3*BL+2*B*D 
TOP = ———___———_—. - -_——__ *B0T 
DSC#¥4*D*(D“+2*B*D) D<+2*B*D 


or making 


A = HM*(B+D)/(DSC#4*D* (D* 


2 


+2*B*D) ) 
AB = (3*B“+2*B*D)/(D“+2*B*D) 
we have 

TOP = A-AB*BOT 
meplacing here 

BOT = AC*TOP 
we obtain 

A 

ToP =-=-__-—_— 
1+AB*AC 
Having already found TOP and BOT we may use Bryansformula to 
determine S;since the thickness of plating is for this case 
smaller at the bottom we must consider this thickness and 


along with it DSC because this will be the larger compressive 


Stress the bottom will be subjected to. 
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2 
k BOT 


From vs Q = * 2 
0, 30396*(1-POI*POL) S 





We obtain 
1 
eens 


1-POI*POT 
0. 30396% —_—————$——* DSC 
k 


or making 





1 
7 = 
L=—FPOLarerL 
0 FOS C6 == er iit 
E 
we have 
So Acre Ouk 


A computer program was written and applied to two parti- 
cular sets of dimensions and two different materials namely steel 
end aluminum, The results are presented next in the computer 
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HOGGING MOMENT AND SAGGING MOMENT OF EQUAL MAGNITUDE 
BiolGN STRESS IN TENSION AND DESIGN STRESS IN COMPRESSION OF 


EQUAL VALUE; LATERAL LOAD ALSO CONSIDERED 


The lateral load being here considered will be a "head" 
of salt water of eight feet on the top plating (HTOP=8 ft) 
and a "head" of salt water equal to the draft on the bottom 
plating (HBOT=DRA). | 

We will use here the simplified approach as suggested 
in the following note; this simplified approach is at least 
conservative, 

NOLeG. blelechy1rom Nis experiments and Cneoretical studies 
ma the effect of normal pressure in the buckling strength of 
thin panels, concluded that the buckling strength is increased 
when combined with a sufficiently large normal load, and that 
the increase is greater for thesimply supported edges situation 
than with fixed edges, For heavier plating , where thedeflection 
does not exceed half of the plating thickness, the increase in 
buckling strength resulting from normal load is negligible, 
me was therefore suggested, for simplicity, that Ene 2errcical: 
stress be computed for in plane loading only, 

The tertiary stresses we are interested in considering 
here are those exerted at points B directed longitudinally; these 
will combine with the longitudinal bending stresses or primary 
stresses as described in a previous ae 


Since now we have DST=DSC and HM=SM=BM we may think of 
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a diagram for the stress distributions as indicated next figure: 


SU souk elie Sulgl Soleas ih 
—> a 


PUP goles): Bh SBS SMEG Gis 


a ne ree 


COMPRESSION TENSION 


From this diagram we may understand that the neutral axis 
should be located closer to the bottom which means that the thi- 
ckness at the bottom will be larger than at the top BOTD>TOF , 
So, we may already write the following formulas where TS was 


replaced by BOT 


DITOP = D*(B+D)*BOT/P 


DIBOT = D*(B*TOP+D*BOT)/P 

I = 4*D*D*((3*B“+2*B*D) *TOP+ (D2+2*B*D) *BOT ) *BOT/P 
where 

P = 2*D*BOT+B*(TOP+BOT) 





see, 


We may now think of the following relationships to deter- 


mine TOP and BOT, 


Sahay 





SBH 
Sioa! = obn-OLTOP 


or knowing that 








STR3T = 0,152222*HTOP* (——) © 
TOP 
STR3B = 0,152222*HBOT* (——)* 
BOT 
and anticipating 
S = Z*TOP 
where 
1-POI*POI 
-Z = 1/( \}0. 30396* —————"*DST ) 
E 
we have 
S 
Z° = (—)? 


TOP 
S 
If we now make = Ze 
BOT 
in order to simplify the equations, we see that when we do this, 


since 
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S 5 S 5 
=a) SS) 
BOT ee 


we are increasing the value of STR3B and this means that since 
DSi — obn+STR 3B 
we are decreasing"slightly" the value of SBH that could be allo- 


wed to reach a larger value; so , the substitution 


S 
(——~)* = 2° 


BOT 
is on the safe side and has the advantage of simplifying the 
equations considerably. 

So we may write 

STH = DST-STR3T 

SBH = DST-STR3B 


Zz 


SER oO. elo 2a22 Rh TOPR*Z 


STR3B = 0.152222*HBOT*Z* 


A = DST-STR3T 
AB = DST-STR3B 
SO 
A*DIBOT = AB*DITOP 
and we obtain replacing DITOP and DIBOT by their equations, 
and simplifying 
AB*(B+D)-A*D 


TOP =————————- * BOT 
Ao 


or making 
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AC = (AB*(B+D)-A*D)/(A*B) 
we finally obtain 


TOP = AC*BOT 


Now In@=Greaer to stress tne material to its design stress 
emtray sClcner COnslder Sli, oo, SEH, SBo. 

Due to ease of calculation - may either choose STS or 
ues oie expressions will be exactly the same, I will consider 


here 


Sukie BM*DITOP/I 


and 


ods DST-STR3T = A 


Substituting for DITOP and I we obtain 


BM*D*(B+D)*BOT 


he D2*BOT* ( (3*B<+2*B*D) *TOP+(D“+2*B*D) *BOT) 


and after simplifying we arrive at 


BM*(B+D) 3*B<+2#B¥D 
BOT = ——————_,»———— - + * T0P 
A*4#D#(D~+2*B*D)  D<+2*B*D 


Making 


AD = BM*(B+D)/(A*4*D*(D<+2*B*D) ) 


2 
AE = (3*B~+2*B*D)/(D“+2*B*D) 


we obtain 


BOT = AD-AE*TOP 
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substituting 


TOE AC*BOT 


we obtain 


AD 
BOT =——_——_—_ 
1+AC*AE 


Considering now the instability formula (BRYAN) we notice 
that since the top will have the thinner plating we will consi- 
der TOP stressed to the maximum allowed DST. 


So we start with 


1-POI*POI 
1/(\10, 30396* ——-———*DsT  ) 
E 


Sie ea) 


N 
I 


(Cp) 
i 


The total stresses will then be given by 


TSTH = BM*DITOP/I+STR3T 
TSTS = BM*DITOP/I+STR3T 
TSBH = BM*DIBOT/I+0.152222*HBOT#S~/BOT@ 
Toeea= BM*DIBOT/I+0, 152222*HBOT*S“/BOr@ 


The computer printouts for the same cases analysed in theprevious 


section are presented next, 
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STRENGHT OF SHIPS#ON THE NEUTRAL AXIS LOCATION PAGE 


1¢ 


30 


4 


5@ 


*AL* IS THE LENGHT 

AL=S00-6 

REAL /Re 

*ORA* IS THE nRAFT 

ODRA=B/3.« 

DA 7 2s 

*«POT* TS THE eEOISSON'S RATIO OF THE MATERIAL 
FOR STEEL 

POTsxAe3 

«oe IS THE YOUNG'S mOpnULUS OF THE MATERIAL 
FOR STEEL 

F=2Dex1 2 een 

*OST* ITS THE nESIGN STRESS 

DST=2@50¢C2Z6 

*PMe TS THE BENDING MaAMENT 
R“54e37143eBuengawAl wAy 

W=0026396*(1+eD0T *POl)*DST/E 

U=SQRT(W) 

Liege, 

eHTOPF* ITS THE HEAD OF WATER ON DECK 

HTOP=HR.e 

*HBOT* IS TRE HEAN OF WATER ON BOTTOM 
HEOT=DRA 

STR3OTHCZe1522 la eH TOP a Zu Z 
STRBEB=HS%*1H82eC ele ae4ROTae lu Z 

A=9ST&STR3T 

ABeNSTeSTR3B 

AC=(AB¥(RtNeguD) /( AMR) 
ADEPM*(B4D)/ (pe ho tDe(N*D420%BaN ) } 

AE =(3e*B we Bele au RQ¥D) JS (Da Dt 2 oe kBHD ) 

*TOP* TS THE rHICKNESS OF THE TOP PLATING 
*80T*® IS THE YHICKNESS UF THE RUOTTOM PLATING 
BOTHAD/(1+4ACwAE) 

TOPREAC#SAOT 

WRITE (5940)TOp 

FORMAT('1tse¢x,'TOPs'F1964) 

WRITE(S220) 607 

FORMAT(! ',2@y,'BOTs'c1@ed) 

*S* ITS THE SPacING OF LONGITUDINAL STIFFENING MEMBERS 
SEZ" TOP 

WRITE(5230@)S 

FORMAT(! 't,2?Px,'S=tFiges) 

Pe2.*D*ehOT+Be( TOP+BROT ) 

*OYTOP* IS THe NISTANCFE OF THE TOP FROM THE NEUTRAL AXIS 
DITOPED* (BaBO7s0"ROT) /P - 
WRITE(5,4A)0170P 

FORMAT(' t,iRy, (DI TOPa2!F4G.4) 

sOTBOT* IS THE DOISTaNCE OF THE BOTTOM FROM THE NEUTRAL AXIS 
DIBOT&O*(BaeTOPsDeROT) /P 

WRITE (5s59)0InQT 

FORMAT (! ',L1Ry,'DIROTatF {Ges ) 

eGMI« IS THE MONENT OF INERTIA OF THE X SECTIOK ABOUT Neds 
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STRENGHT OF SHIPSRON THE NEUTRAL AXIS LOCATION PAGE 2 
GMT =heeDeDe ( (9, *#8*R 42, *B aD) © TOO *BOTH (0*D042.*BwL ) eROT* BT) sp 
‘@ #TSTH*®IS THE eTRESS ON TOP DUF TO HOGGING 
e *TSTS*IS THE cTRESS ON TOP DUS TO SAGGING 
© «TSBH*IS THE stTOeSs On BOT DUE TO HOGGING 
Cc xTSBS*IS TRE «erRESS Oy BOT DUE TO SAGGING 
C ¥#STCRIT* YS Tur CRITICAL STRESS ON TOP 
‘@ *SBCRIT* IS Tuk CFITICAL STRESS ON BOT 
TSTHERM*DITOP ,aMresrRaT 
TSTSERBM*DITOP/GMI4STRAT 
STCPIT=TOP*TOn#F/(4.37396*S*Se(LewPOIT*PO!] )) 
WRITE (5260)TSTuH | 
Ban FORMAT(' 'y19y,!'TSTHEIF1Ge3) 
WRITE (52 70)TS7S 
7a FORMAT(' ',29x,°TSTSEVF1G03) 
WRITE(S280)STCRIT 
Ba FORMAT(' 's17y,'STCRITE'FIF®3) 
TSRBHEBMXDIBOT/GMI+2 1522202 eHB OT *¥S¥*o/(BOT*BOT ) 
TSBSEBM*NIBOT/GMI th el RP2Q2KHROT*SHSS(BOT*BCT ) 
SBCRITEROTHROTHE/S (A. 3G39648*S a (LeePOlePOl]) ) 
WRITE(S*90)TSRH 
9A FORMAT(' ',19y,!TSRNEIF IA 2) 
WRITE(52100)TSAS 
17G FORMAT(' 's19yv,'TSAS=HF1 Ge?) 
WRITE(S#11A)SRCRIT 
41@ FORMAT(' 'y17yv,'SBCEI rare 63) 
END 
(S) eV @Eé4CCV) AL 2654 lv) B A6RCCIV) DRA A664CVI OD 
Oty) Por A674CV) E AAAACS)] eR C6R84CV) DST ZS8CCV" BM 
BetV) WwW QFAQCVI VY ABAAlLS) SQRT A6BA4CVI Z ABARBCV* HTOP 
may) HaoT ASBACV) STR3T AHRRCY) STRIB C6RCCVI A Z6CO@CV! AB 
V) Ac @éDR CV) AD ABECCy)] AE C6FATV] BOT 26F4CV2 TOP 
tfL) 12 GAPACS) al PA3aT] 20 C6FRCVI S QA7CEL® 3a 
(Vv) P Q7A9@CVI OTOP Pees lr) 42 @7a4CV) DIROT A3Z54CL9 57 
RCV) Gm] a714CV2 TSTH A7iR Cy) TSTS G@7iCCV) STCRIT @464CL°: 6f 
{L3 7o P4NACLI RO 77eaCy] TSBH C724CLV)] TSas A72RCVI SRCRIT 
atL) 997 GOFBEELJ 100 26ect J] 11¢ CAGACS] eV 
TOGQAM | ABELS! 
2O7A suaAtIne 2ERe eV 2CEe ecOMP 2984 ALOG 2ECe @] (ag be OE 
2DAA e7ERO CD7R $6 270C a4 2E26 .MES 2CBH oW 2A8C E 
2Ci2 ANT 28BE SORT 2D4A #RARG 2NA& 45 eAsC AEXP 2DAE 
@07C $g 2DAE +O 2E84 oy 3EDC 
ITRY@POYNTs: 
F73C eR C7DC oA PRRE SART 2984 AaLOG 2ARC AEXP ePaSC € 
eCi2 AYNT CCR ow 2CcEr ecOMP PHS £6 CD4A eRARG eo7c $ 
2046 eS 2DAA «ZERO POAE ef RCNT erBba .0 CEBE MES PEIA 
2EB6 ey CECA al 





F 


_— 
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BOT= @.8283 
S= 11.,a@49 
DITOP= 35,9820 


DIBOT«= 20.4736 
TSTH= 2509%¢0.%1?P 
STS SP5hCAA.A12 
STCRIT= 2@5A°%Aa HG 
TSBH= 1147P2.468 
TSBS= 147294458 
SBCRIT=7°66365,2524 


IO OPERATING SYSTEM VEReETION 1 REVISION Ale AI/A4/74 GENERATED (4727/74 


HANDLING CHARGE $ «35 / JOB ¢35 
LINES PRINTER PR1 S$ 1475 / K LN °14 
CARDS READ S lega / K Ch °13 
PENTTeE? VECTORS S$ 25 / {ARF oA? 
MONEL 77 SECONDS $25-¢aq / HOUR 049 
MOREL kx SECONDS S420eRnm 7 HOLE 0An 

TOTAL CHARGE & 074 


49% 44731 LOGGED OUT g4729/74 18:3:35+* % 33°19 LEFT AFTER 42 LOGINS: 
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Note here that because we have in this run BOT>D> TOP the 


approximation made to simplify the formulas, namely 
5 
7, = 





ekOME 
becomes very inaccurate compared to what the actual value 


should be 
3 





yp 
Or 


this is the cause for obtaining such low stresses at the bottom, 
We might refine the result for this case as is done for the 
"1000 ft steel case" also included in this section which presen- 
ted still another difficulty, Since all steps to have the bottom 
stressed to DST are explained in detail for such case they will 


be omited here, 
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E@eomedsceecavcwamresc il CllOse LO Lhe correct solution with 
mac Hovtom stressed to almost DST, 

Next is presented the "1000 ft steel case" where still 
another difficulty was encountered and a procedure to have the 


bottom stressed to DST is exemplified, 
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STRENGHT OF SHIPS#ON THE NEUTRAL AXIS LOCATION PGi 





10 


2B 


30 


40 


52 





#AL* IS TRE  -ENGHT 

ALH=1l@Cve 

PFAL/5Se75 

wDRFA® JIS THE DRAFT 

DRA=H=B/3°3 

D=AL/14>° 

*FOI« IS TRE POISSON!S RATIO OF THE MATERIAL 
FOR STEEL 

FOT=HC+3 

Beers THE YOUNG'S MODULUS OF THE MATERIAL 
POR STEEL 

er yPem iia ® +6 

eC ST* IS THE DESIGN STRESS 

DETHSPOV Ke 

*Em* IS THE BENDING MOMENT 
Bk=1e37143*eR*ODRA»AL#AL 

w=Ge 3U396"*(12.9°P01*«PGI)*OST/E 

U=SORT(h) 

Z=1e/U 

emTOrP=s [S THE HEAD OF WATER ON DECK 
HTOP=Ke 

*hROT* IS THE HEAD OF WATER ON BOTTC 
WEOTSDRA 

ei Rete «1 See ee *HIiCP=Z+Z 
STR3BB=HE4 61 SPEPAO#*HROTe Za Z 

SapST=SIRST 

AE=DbST-STK3B 

ACH=(AB4*e(BR+C )eAx*D)Y (ARR) 
AD=BN*(B4+D)/(Ae4eete (D*N+26*Hl)) ) 
AEH(BSethe b+ lek aN) S(DaeD+l ou hsD) 

*TOPs IS THE THICKNESS OF THE TOP PLATING 
*BOT* 3S THE THICKNESS OF THE BOTTOM PLaTING 
BOT=ACD/(1¢e+AC*AE) 

TOP=AC*8OT 

Weerlet 5, 1) TOP 

FORMAT('1',evXs'TOPH'F1P 64) 

WRITE (S220) 86T 

FORMAT(' ',27X,'BOTH'FIC +4) 

Poeelomine oe ACTING OF LONGITUDINAL STIFFENING MEMBERS 


S=Z*TOP 
Wreewe(5,30)S 
FORMAT(' tyerpx,'Sa='F1H¢4)} 


P=Pe*DeBCT+H a ( TOP+E CT) 

*DITOP*« [S THE DISTANCE OF THE TOP FROM THE NEUTRAL AXIS 
DITOP=0*(6«BOT+D«e ROT) SP 

BRT TE(5246)DITOP 

FORMAT(! ',138X%,'OTTOPH=IF IAS) 

*“OTBOT* ITS THE DISTANCE OF THE BOTTOM FROM THE NEUTRAL AXIS 
SIBNT=N*(BeTOP+0"50T) /P 

Vink iere| SHES wDIBOT 

FORMAT(' ',18X,'DIBOTH'Fide4d) 

*GMI*® IS THE MOMENT OF INERTIA OF THE x SECTION aBUOUT Nehe 
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STRENGHT OF SHIFS*®ON THE NEUTRAL AXIS LOCATION PAGE ? 
2 GMT E4 ot OuDe( (36 *Aw A+ Pe KAHN) *TOP*ROT+ (Dali +2 ee BD) mBOT*ROT)/P 
E¢ Ee «TSTH#*1S THE STRESS ON TOP DOUE TO HOGGING 
C *TSTS*¥IS Tee STRESS ON TOP DUE TO SAGGIAG 
C »TSOH*® {TS TRE STRESS ON ROT DUE TO HOGGYAG 
a .'obSe7o (=e ommeosS ON ROT OWE TO SAGGING 
C SOPCRIT*® IS THewGRITICAL STRESS ON TOP 
C ~SRGRIlTT4# IS.ime GRITIGAL STRESS ON KGT 
TSTHEBM*DITOPSGMI4S7TR AT 
TSTS=AM#*DITOFSGMI+S7 RAT 
STCRITHETOP¥TOP*E/( 2.30 396"Se5e(1e"POT*POT) ) 
g WRITE (Ss 6&6) TSTro 
b A OG FORMAT(' ',1lSxXj 7 SG 1eatr id «- ) 
TE wWRITE(S*#70)TSTS : 
ST 7 PGRmAT(! “5 Poxe 8ISTS=tt10e 3) 
5 We i LE | Gore tl) Sea 
D BO FORMAT(' ',17X,*SICRIT='F10¢3) 
TSBHERMS DCI ROT/GNI4+¢.1 5°20 22sHKHOT*S*&S/(BOT*ROT ) 
TSBSH=AM*FOLSOT/GMIT 47.4522 02C0*¥HROT*Sa*S/(BGOT*B0T ) 
SECRITHBOL+*BOT*Es( 4.32396 "S*Sue(LemPOL*Pol) ) 
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It does not come aS a Surprise that the program does not 
work for the "1000 ft steel case" since in this situation we 
obtain STR3B DST which cannot be accepted and also this 


program is not prepared to handle, 


AL = 1000 
B= AL/5.75 = 173.913 
DRA = B/3.3 = 52.7 
POk-=7 0.3 

E = 30*10° 
Dote—e2000 


BM = 1, 37143*B*DRA*AL“ = 12569,45*10° 


W = 0, 30396*(1-POI2)*DST/E = 230.503*107° 
15.182#107? 
S5), oy 


HBOT = DRA = 52.7 


U = W 


Z = 1/U 


STR3B = 0,152222*HBOT*Z~ = 34803.57 DST 
Since we understand that the program should work for DST>STR3B 
we make the next run with DST = 35000 to prove that in such case 
a solution is obtained,and after proceed to obtain the actual 


solution for the case we were considering with DST = 25000, 
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Now let's consider our situation with DST = 25000. Once 


the situation is such that due to the approximation made in 
S 


————— = 7 
Boe 
we end up with a STR3B> DST we must consider a value Z1 lower 


than Z to bring STR3B down, This can be also understood since 


it is acceptable to make here 
S 
ee = =, 


TOP 
but because we obtain with the program as is TOP<«BOT, making 


the approximation 

mune = /, 

BOT 
makes Z attain a much larger value than it actually should, 
‘We must then consider Z1<Z = 65.8 and may start by trying 
Z1 = 35 


The computer printout in this case is as follows, 
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This time we obtained a realistic but still unacceptable 
result; we notice that the bottom is being stressed to 31000 psi 


when it can only go up to 25000 psi; besides we assumed 








Zi = = 35 
BOT 
and end up with 
= 225.4 
Se = NY UG 
BOT on OM 


This means that the Zi value should be higher, 


We may now try 421 = 40; the computer results for this case are 


shown next, 
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In this case we assumed Z1 = 40 and ended up with 





S Lo) 6 D 
eee TT 7 
BOT 6.59 


and a stress at the bottom of only 19000 psi which is too small, 
We may now try to use a better Z1 value and can do the 


following: 
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Stress at the bottom 


5*6024 
Zl = 35+ = Sh? (eats 
Avono 
20000 a 
Le 
41 ASSUMED 


4o 
'we will use this Z1 value 
, The results obtained with this Z1 value are the following: 
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The results show that we have obtained 


S AN By ian 


while it was assumed 

tN = 3/017 
We could further correct the result until Z1 assumed and 241 actual 
mould coincide but 10 seems Chat there is no need to refine more 


the solution for the purpose here involved, 


Note that this procedure should also be followed for the other 
cases in order to bring the stress at the bottom to its maximum 


allowed DST for such cases as well. 
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This 1000 ft aluminum case shows again a stress at the 
bottom much lower than DST, It should be corrected as done for 


sev0O it Steel” in the previous Situation. 
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HOGGING MOMENT LARGER THAN SAGGING MOMENT 
Piolo SIRESS IN COMPRESSION LARCER THAN DESIGN STRESS IN TENSION 


LATERAL LOAD ALSO CONSIDERED 


With the help of what was understood for the similar si- 
tuation without the lateral load being included, we may here 
proceed at once to look at what the diagram of stress distribu- 
tions might look like. Actually, the diagram that will be con- 
sidered first is basically the same with only a small change 
Gue to the addition of the tertiary stresses; since these will 
probably be of the same magnitude on both top or bottom and 
assuming that their actual value does not differ considerably, 


what we may start looking into is a diagram as follows: 


DSC | DST 
Cz 
| : 
| | \ oa) 
| ; 
; | | Wz ; | 
ee \Le bee’ 
pie =: KE aan ee as 
| | va I 
| 
, VA , | 
| a ' i 
Sere | <2 
stat —S'7R3 


We would consider here the same basic distribution of longitu- 


dinal bending stresses thus obtaining 
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STS STH STRBL 
ie { 
| Ss! 
| ue | 
a : DEEP 
: Ea | 
{ JL NA. 
en pes Sg i I A. 
I Vee i 
¢ 
| / ) 
1 Y | DIBOT 
Va ; , 
| wv 
-— 
STR3B SBA SBS STR3B 
; D or 


We Will then expect to obtain TOP> BOT and DITOF< DIBOT 
and we may already consider the respective expressions for 
DITOP, DIBOT and I with TS replaced by TOP 

From the diagram we may consider thefollowing equations: 

Slo reOr=SBH=DITOP 
and notice here that if we want to stress the material to the 
allowable design stresses we have 


STH 


DST-STR3T 


SBH DSC-STR3B 
we may combine these equations to obtain 
(DST-STR3T)*DIBOT = (DSC-STR3B)*DITOP 


STR3T and STR3B are given by 


S 
STR3T = 0,152222*HTOP*(——) ~ 
TOP 
S 2 
STR3B = 0,152222*HBOT*(——~ 


BOT 





anGrcrpaving here that S = Z*BOT because BOT is the thinner 


plating we have 


and thus may write 


STR3B = 0,152222*HBOT*Z* 
Since 
S S 
(ere (ye 
TOP BOT 


because TOP> BOT, if we make 


3} 
4a 
TOP 


Zz 


to simplify the equations that will follow we will understand 
that in doing so we are only increasing the value of STR3T, 
This as we see fron 

STH = DST~STR3T 
will make us have STH with a slightly lower value than the one 


that is allowed, so we are on the safe side when we consider 


a 

TOP 
Now if we end up too much on the safe side and we have a TSTH 
too low we may bring it to a value closer to DST by considering 
@Z1 value as done and exemplified in the previous section and 
done again in this section for the case 1000 ft steel where 


for this particular exemple it is TSBH that needs correction 


as explained in detail 


ll 
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We may then write 


1 


1-~POI*POI 
0, 30396*%*————————-*DSC 
E 


ie 


a 


STR3P = 0,152222*HTOP*Z 
STR3B = 0,152222*HBOT*Z~ 
and from 
(DST-STR3T)*DIBOT = (DSC-STR3B)*DITOP 
we have if we make 
RS SPSS INS Shh 
AB = DSC-STR3B 


P = (2*D+B)*TOP+B*BOT 


D*TOP*(B+D) D* (B¥BOT+D*TOP) 
A ce 
le ie 


which gives after simplification 

A*{(B+D)-AB*D)*TOP = AB*B*BOT 
and making 

AC = (A*(B+D)-AB*D)/(AB*B) 
we have 

BOT = AC*TOP 

Following th e suggestions of a previous section we will 
also use the equation that will make the bottom to be stressed 
the most in compression under the larger of the moments being 
considered (HM in this case), 
We obtain 

SBH = HM*DIBOT/I = AB 


Replacing DIBOT and I by their expressions we obtain 
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HM*D*TOP*(B+D) 
Sf — i a Se ee 
Lap <*7OP*( (3*B“+2*B*D ) *BOT+(D“+2*B*D) *TOP) 


simplifying we end up with 
HM*(B+D) 3*B“+2*B*D 
Tor =+—__-——___—_—  - ae Oe 
AB*4#D*(D<+2*B*¥D)  D<+2*B*D 
now making 


AD = HM*(B+D) /(AB*4*D* (D“+2*B¥*D) ) 


AE = (3*B°+2*B*D)/(D°+2*B*D) 
we have 

TOP = AD-AE*BOT 
substituting here 

BOT = AC*TOP 
we obtain finally 

AD 

TOP =-———__—— 
1+AC*AE 
Considering here again the bottom as the location of the thinner 
plating (in this particular situation we are looking at), and 
knowing that it will be subject to a compressive stress of 
magnitude as large as DSC we may again write S$ = Z*BOT 
We finally obtain the total stress on top and bottom in both 
tension and compression by adding the longitudinal bending 
stresses to the tertiary stresses, 


TSTH HM*DITOP/I+0, 152222*HTOP*S“/TOP* 


Arse s 


SM*DITOP/T+0. 152222*HTOP*S</TOP* 


TSBH = HM*DIBOT/I+STR23B 
TSBS 


SM*DIBOT/I+STR3B 
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NOs SouticecmUslally ~hewsceonuary Stresses are assigned 
a constant value like 2000 psi or 3000 psi they can be a- 


ecounted Tor when giving a value vo the design stresses, 
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AS We See the results obtained for the 500 ft aluminum 
case are acceptable and agree with all our considerations in 
the derivation of the formulas, namely TOP> BOT, 

The foliowing two runs, for 500 ft steel and 1000 ft alu- 
minum show that the results are in disagreement with TOP> BOT 
used in the formulas, We will see next why this happens and 


derive the proper set of formulas to be used for such cases, 
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STEEL 
role= "Ors 

E = 30*10° 

AL = 500 

B = AL/& = 62.5 
DRA = B/3 = 20,833 
D = AL/9=55. 556 


1. 37143*B*DRA*AL* = 44642,187*10 
Ly 


HM : 


SM = B¥DRA*AL“ = 32551. 562*10 


W = 0, 30396*(1-POI“~)*DSC/E = 276. 6*107° 

Z = 60,129 

eRe 6 

Heel = 205633 

STR3T = 0. 152222*HTOP#Z” = 4402, 865 note that this value 
should be corrected for the 4 
approximation and actually is 
Slightly smaller 

STR3B = 0,152222*HBOT*Z* = 11465,61 


We see from the tertiary stresses obtained that we end up with 


the following diagram of stress distributions: 


20597 —alt9n. 
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From the diagram we see that theneutral axis will end up closer 


to the bottom and so we must consider BOT> TOP in the formulas, 


The same type of diagram will be obtained for 1000'aluminum case 
Since the tertiary stresses obtained in Such case are 

1498, 73 

9873078 


Sno 


Sle Sie. 


For the 1000! steel case the same type of diagram will occur 
due to the tertiary stresses being 


STR3T = 4402, 86 


il 


ik 


Sos 29004, 

Besides this situation requiring us also to consider the 
neutral axis closer to the bottom, still another problem arises 
since STR3B is larger than DST=25000, This case is solved by 


properly correcting the inaccuracy that was introduced by the 


Z approximation as exemplified next. 
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Now let's derive the formulas capable of handling the 500! 
steel case, the 1000! steel case and the 1000' aluminum case, 
We must consider BOTD TOP and DITOP> DIBOT 


DITOP = D*(B+D)*BOT/P 


|OUL leiOMk D* (B*TOP+D*BOT ) /P 


I = 4*p2*Bot+( (3*B“+2*B*D)*TOP+ (D2+2*B*D)*BOT)/P 


Wwnoere 
P = 2*D*BOT+B*(TOP+BOT ) 
Starting again with the equation 
Ske DIBOT = seo *DITOP 
and making 
A = DST-STR3T 
AB = DSC-STR3B 
we obtain 
A*DIBOT = AB*DITOP 
Knowing that now the top is the thinner plating we anticipate 
Se= 2 er 
and may write 
STR3T = 0,152222*HTOP*Z“ 


Zs 


Sulye Sis See 22 NeOlLaZ 


where the Z value in STR3B should be which has a smaller 





BOT 


value than Z, 
We then obtain: 


A*D*(B*TOP+D*BOT)/P = AB*D*(B+D)*BOT/P 
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Simplifying and making 
AC = (AB*(B+D)-~A*D)/(A*B) 

we obtain 
TOP = AC*BOT 

Making now the top to be stressed the most in tension we have, 
STH = HM*DITOP/I =A 


Replacing DITOP and I by their equations we have 


HM*D* (B+D )*BOT 


4eD<*BOT* ( (3*BO+2*B*D) *TOP+(D2+2*B*D) *BOT) 


and obtain 


HM*(B+D) 3%B<+2*B¥*D 
ee ee i ee ee eT OP 
A¥4%D* (D£+2*B*D) D©+2*B*D 
Making 
AD = HM*(B+D) /(A*4*D* (D©+2*B*D) ) 


AE = (3*B2+2*B*D)/(D“+2*B*D) 


the formulas become 
BOT = AD-AE*TOP 
and substituting 
TOP = AC*BOT 
we finally have 
AD 
BOT = 
1+AC*AE 


The results obtained for the examples chosen are presented next: 
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We could have obtained TSBH = 30000 by correcting the 
approximation used Z = S/BOT as is done for the 1000! steel 
case that follows next, 

In this case we already anticipate the use of a Zl 
value as done in the previous section and according to the 


output obtained we will then obtain the correct result. 
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Now we see that with an assumed value Zl = 35 we obtain 

S = 215.1808 

BOD =" 6.60 so Z1 actual = 58, 38 

TSBH = 47512.793 (too high) 

ELSES also too high 

The Z1 actual value obtained suggests that we should have 
assumed an higher Z1 value and with an assumed value Z1 = 40 
we then obtain 

S = 203.626 

BOT = 4,428 so Z1 actual = 45.98 

Toph eot lol cog | voolhiah) 

Lope also too high 
This suggests an assumed Z1 value higher than 40 but now we may 
interpolate with the two results already obtained and determine 
the Z1 value that will bring TSBH = DSC = 30000 psi and done 
properly this new 41 assumed will end up being Z1 actual, 
From interpolation we obtain 

(40-35) *( 34131, 855-30000) 
Zt. = 404 ————_—— = 41, 44 
(47512, 793-34131, 855) 


and now assume this value obtaining the following result: 
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We see that TSBH is just a little too high and the same 
Happci omar toes, Mere 1t mignt be the case that TSBS should 
have been the one under consideration for the choice of a Zl 
value. Nevertheless for the purpose of this work this example 
. Snows the procedure to be used along with the equation derived, 
The accuracy obtained can be total. Here we assumed 
Z1 = 41,544 and ended up with 41 actual = 41,83 which is close 
enough for our purpose and TSBH and TSBS is very close to what 


it should be to use the material effectively. 
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This last example shows that TSBH is low and the Z 
approximation used in the formulas here again should be 
corrected as already exemplified in order to bring it closer 


to the DSC value, 





234 


CONCLUSION 


Heenecnis Work the author ws led to conclude that there 
is no simplified approach to the neutral axis location in the 
midship section, capable of handling properly all or only 
a majority of cases, It is advisable to look into each parti- 
cular situation by first gaining some insight at what the 
diagrams of stress distributions may look like as exemplified 
in several sections, The two last sections exemplify well why 
this may be of some advantage and specially the last section which 
shows that while at first it was expected for the examples stu- 
died a neutral axis closer to the top, actually the results 
ended up requiring a neutral axis closer to the bottom and a 
new set of formulas had to be derived, 

As a future Aeveveenent for this work it is here sugges- 
ted that a more ship like cross section could be analysed as 
a step forward towards obtaining a method of determination of 
scantlings as direct as possible, This next step should consi- 
lone Ocrivatien ot proper formulas for DETOP, DIBOT and I 
involving only the plate thicknesses as unknowns, It should 
be considered then, in order to simplify the problem., a stan- 
dard type of frame and its scantlings should be constant throu- 


gehout the midship section, 





J.H, EVANS, 
ie EVANS, 
J.H, EVANS, 
J.H, EVANS, 
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